CYCLICALLY COUPLED SPREADING AND PAIR ANNIHILATION 



HAYE HINRICHSEN 



Abstract. Recently it has been shown that the transition of the 1+1-dimensional annihilation-fission 
process 2X— >3X, 2X— >0 exhibits an unusual type of noncquilibrium critical behavior. The phenomeno- 
logical properties of critical clusters are characterized by two different dynamic modes for spreading 
and diffusion. In order to describe the interplay of these modes, we introduce an effective model which 
involves two species of particles A and B. The A-particles perform an ordinary directed percolation 
process while the B particles diffuse and annihilate. Both subsystems arc cyclically coupled by particle 
transmutation A <-> B. The resulting critical behavior is in many respects similar to the one observed 
in the annihilation-fission process. 



1. Introduction 

The study of phase transitions into absorbing states is a fascinating field of nonequilibrium statistical 
physics 0] . Such phase transitions can be observed in models for the spreading of some (generally non- 
conserved) agent as a result of a competition between local reproduction and decay. If the rate for 
reproduction is sufficiently high, the system is able to maintain a fluctuating active phase where the 
stationary concentration of the spreading agent is positive. On the other hand, if the decay process 
dominates, the concentration of the spreading agent decreases and tends to zero. Eventually the system 
becomes trapped in an absorbing state from where it cannot escape. An interesting situation emerges 
at the borderline between survival and extinction. Here the system undergoes a nonequilibrium phase 
transition which is characterized by non-trivial critical behavior. It is believed that transitions into 
absorbing states can be categorized into a finite number of universality classes. Typically each class is 
associated with certain symmetry properties of the dynamics. 

At present two universality classes are firmly established. The first and most prominent one is the 
universality class of directed percolation (DP) The DP class is characterized by the absence of 

symmetries (apart from conventional symmetries such as translation and reflection invariance) and covers 
a wide range of models d,||. Roughly speaking, DP models follow the reaction-diffusion scheme X — > 2X, 
X — > 0. In addition, there has to be a nonlinear mechanism which limits the particle density. In 
'fermionic' models with at most one particle per site, the density is limited automatically. In 'bosonic' 
models, allowing for infinitely many particles per site, this mechanism has to be implemented by adding 
the reaction 2X — > X. It is important to note that DP is a unary spreading process, i.e., individual 
particles are able to reproduce and destruct themselves. 

Some time ago Grassberger et al. discovered a second universality class which can be considered as 
a generalization of directed percolation from one to two absorbing states related by an exact ^-symmetry 
(DP2). It comprises various models, including nonequilibrium Ising models ||, certain monomer-dimer 
models |kJ , as well as generalized versions of the Domany-Kinzel model and the contact process [[ll] . In 
1+1 dimensions it is possible to regard kinks between differently oriented absorbing domains as particles. 
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By definition, these particles evolve according to a parity-conserving dynamics, exemplified by branching- 
annihilating random walks with two offspring X — > 3X, 2X — > ]T^-[r^|. For this reason, the universality 
class is often referred to as the parity- conserving (PC) class. However, it should be noted that the PC 
class and the DP 2 class are different in higher dimensions. 

Three years ago, Howard and Tauber [jl5| raised the question whether there might exist a third 
universality class of models with a phase transition from 'real' to 'imaginary' noise. As a prototype, 
they introduced the annihilation- fission (AF) process 2X — > 3X, 2X — ► with single-particle diffusion. 
Obviously, this process is neither parity-conserving nor invariant under any other unconventional sym- 
metry transformation. In contrast to DP and DP2, the AF model is a binary process, i.e., only pairs 
of particles can decay or reproduce themselves. The critical properties at the transition are still poorly 
understood. Recently, Carlon et. al. analyzed the transition using density matrix renormalization group 
techniques Jig] . Estimating the critical exponents they arrived at the conclusion that the transition of 
the AF process should belong to the DP2 universality class, although there is no Z 2 symmetry or parity 
conservation law. However, as pointed out in Ref. ]l7] , various physical arguments suggest that the AF 
process might belong to an independent universality class which has not been investigated before. 

Since the universality class could not be identified so far, it is of interest to understand the most 
salient features of the transition in the AF process from a descriptive point of view. To this end, we 
introduce an effective model which separates the dynamics of pairs and solitary particles in the AF 
process by introducing two species of particles A and B. The A's perform an ordinary DP process while 
the -B's are subjected to an annihilating random walk. Both subsystems are cyclically coupled by particle 
transmutation. It is shown that this three-state model exhibits a nonequilibrium phase transition which 
is similar to the one observed in the AF process. 

The following Section briefly summarizes the phenomenological properties of the AF process. In Sec. | 
we demonstrate that the critical behavior is governed by two competing dynamic modes for spreading and 
diffusion. Based on this interpretation we introduce an effective model involving two species of particles 
(see Sec. [|). In Sec. || the critical exponents of this model are estimated by Monte Carlo simulations. 
The article ends with several concluding remarks in Sec. ^. 

2. Phenomenological properties of the annihilation-fission process 

Assuming the usual scaling picture for phase transitions into absorbing states, we expect the AF process 
to be characterized by four critical exponents /3, /3', v±, and m. The first one is associated with the field- 
theoretic annihilation operator and describes the behavior of the stationary density p s tat ~ (p~Pc)^ close 
to the transition. The exponent (3' is associated with the creation operator and plays a role whenever 
initial conditions are specified. For example, the survival probability of a cluster grown from a single 
seed involves the exponent 5' = (3' ■ The other two exponents are related to the spatial and temporal 
correlation lengths £j_ ~ \p — p c \~ Vl - and £|| ~ \p — p c | _1/|1 , respectively. It should be noted that in the 
case of DP the two order parameter exponents (3 and /?' coincide because of a duality symmetry (see e.g. 
Ref. |Q). In the present case, however, they turn out to be different. 

Using a density matrix renormalization group approach, Carlon et. al. estimated the critical exponents 
of the AF process for various diffusion rates 0.1 < d < 0.2. Since their estimates z = 1.73 ... 1.81 and 
j3/v± = 0.46 . . . 0.5 were in fair agreement with the numerical values of DP2 exponents (see Table 1), 
they concluded that the AF process should belong to the DP2 universality class, although there is no Z 2 - 
symmetry or parity conservation law in the system. This conclusion, however, collides with the general 
believe that critical phenomena are determined by their symmetry properties or, equivalently, by their 
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class/model 


a 


8 


z 


8' 


V 


DP 


0.2765 


0.1595 


1.581 


0.1595 


0.3137 


DP2/PC 


0.92(2) 


0.286(2) 


1.74(2) 


0.286(2) or 


or 0.286(2) 


AF d=0.1 


0.57... 0.62 


0.25... 0.29 


1.67 ... 1.83 


0.12... 0.15 


w 0.10 


AF d=0.5 


0.38... 0.42 


0.21... 0.23 


w 1.7 


« 0.145 


w 0.23 


AF d=0.9 


0.38... 0.40 


w 0.20 






w 0.48 


present model 


0.38(6) 


0.21(1) 


1.75(10) 


0.15(1) 


0.21(2) 



Table 1. Estimates for the critical exponents in 1+1 dimensions. 



associated field theories. Without the required symmetries on the microscopic level, a special mechanism 
would be necessary in order to restore these symmetries effectively on large scales. As there is no such 
mechanism, it is near at hand to expect that the AF process does not belong to the DP2 class. To support 
this point of view, preliminary Monte-Carlo simulations were presented in Rcf. [fL8| and later improved 
by Grassberger and Odor [l9[] . Because of strong corrections to scaling, the estimates for the critical 
exponents depend on the numerical effort. The estimates for 8 = Vx/vu, for example, are scattered over 
the range 0.25 . . . 0.29 and seem to decrease with increasing simulation time. A similar tendency was 
observed for the exponents z and j3. Moreover, it is not yet fully clear to what extent the exponents 
depend on the diffusion rate d. A tentative list of critical exponents, including recent results obtained by 
simulations on parallel computers [ p0[ , is given in Table 1. 

In order to understand the transition in the annihilation-fission process from a phenomenological point 
of view it is helpful to analyze the spatio-temporal structure of critical clusters. To this end we introduce 
a novel type of scale-invariant space-time plot which can be used to visualize the scaling properties of 
critical clusters in systems with absorbing states. Starting with a localized seed (a pair of particles) 
at the origin and simulating the spreading process up to 10 6 time steps, the rescaled position of the 
particles x/t 1 ^ is plotted against log 10 t, where z = V\\/i>x is the dynamic exponent of the process under 
consideration (see Fig. [l]). By rescaling the spatial coordinate x, the cluster is confined to a strip of finite 
width. Compared to linear space-time plots the scale-invariant representation of a cluster has several 
advantages. On the one hand, it is possible to survey more than four decades in time. On the other 
hand, the plot provides a simple visual check of scaling invariance. Roughly speaking, scaling invariance is 
fulfilled if the cluster's appearance is time- independent, i.e., spatio-temporal patterns should look similar 
in the upper and lower parts of the figure. It is needless to say that this visual check does not replace 
an accurate quantitative analysis. Nevertheless such a scale-invariant plot may improve the intuitive 
understanding of the asymptotic long-time behavior and may also help to identify relevant and irrelevant 
contributions. 

Let us first consider a critical cluster of a branching-annihilating random walk with two offspring (see 
left part of Fig. |l|). Obviously, this process is characterized by an ongoing competition between particle 
reproduction and decay. Contrarily, the annihilation-fission process admits undisturbed random walks 
of solitary particles over long distances. As shown in the middle of Fig. [|, this leads to a very different 
visual appearance of the cluster. It is important to note that these differences persist as time proceeds. 
However, if both processes were to converge to the same type of long-range critical behavior as suggested 
in Jj"6| ], we would expect the clusters to become increasingly similar in the lower part of the figure. As 
there is no indication of such a convergence, Fig. |l| supports the viewpoint of Refs. |l7],[2(]] that the AF 
process might represent a new type of noncquilibrium critical behavior. The corresponding universality 
class should be characterized by the following properties: 
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*~ x/t 1/2 ~< 1 x/t 1/z 1 x/t 




BAW2 annihilation-fission cyclically coupled model 



Figure 1 . Critical clusters generated from a single seed at the origin in a scale-invariant 
representation. The rescaled position x/t 1 /* of the particles is plotted against time on 
a logarithmic scale. The curved lines mark the positions x — ±f0, ±f00, and ±f000, 
respectively. Left: Branching-annihilating random walk with two offspring (BAW2). 
Middle: Annihilation-fission process (also called pair contact process with diffusion). 
Right: Cyclically coupled spreading and annihilation processes as an effective model for 
the AF process. The figure is explained in the text. 
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DP BAW2 AF process 



Figure 2. Cartoon of three different 1+ f -dimensional spreading processes. Solid 
(hollow) circles denote reproduction (decay) of particles. Left: Directed percolation. 
Middle: Parity-conserving branching annihilating random walk with two offspring. 
Right: Schematic drawing of the annihilation-fission process where spreading occurs 
in avalanches (represented as triangles) followed by diffusion of solitary particles (see 
text). 

1. Single particles diffuse but do not react. 

2. Reproduction requires two particles to meet at neighboring sites. 

3. Particles are removed if at least two particles meet at neighboring sites. 

4. There is no unconventional symmetry (such as parity conservation). 

5. There is no frozen disorder in the system. 

6. There is a mechanism limiting the density of particles. 

Consequently, many other processes, such as the fission-coagulation model 2X — » 3X, IX — > X and the 
reaction-diffusion process IX — > 3X, 3X — > X are expected to exhibit the same type of nonequilibrium 
critical behavior. 

3. Interpretation as a spreading process with two species of particles 

To what extent is the transition in the AF model different from ordinary DP and DP2 transitions? As 
Fig. ^ suggests, there are two separate dynamic modes for spreading and diffusion. The spreading mode 
is characterized by sudden avalanches with a high density of particles. Here the dynamic processes are 
dominated by interacting pairs of particles. Once an avalanche has stopped, the system enters the diffusive 
mode, in which solitary particles perform a simple random walk. When two of them meet at neighboring 
sites, they may release a new avalanche, as illustrated in Fig. |^. The asymptotic critical behavior at 
the transition will depend on the relevance of the two dynamic modes. In principle there are three 
possibilities. If the spreading mode becomes dominant we expect a crossover to DP. Conversely, if the 
diffusive mode governs the asymptotic regime, we expect a purely diffusive behavior with the dynamical 
critical exponent z = 2. However, Fig. ^ strongly suggests that both modes are equally important and 
balance one another as time proceeds. In fact, in the rescaled representation the typical spatio-temporal 
patterns do not change over four decades in time. 

In order to investigate this transition in more detail, we suggest a phenomenological explanation of the 
observed spatio-temporal patterns. The basic idea is to describe the two dynamic modes in terms of two 
separate reaction-diffusion processes involving two different species of particles A and B. The spreading 
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Figure 3. Cyclically coupled directed percolation and annihilation. 



mode is governed by the dynamics of A-particles. Roughly speaking, the A-particles can be thought of 
as representing pairs of particles in the original model. Obviously, there is no parity conservation law 
for the number of pairs, hence the ^-particles in the new model evolve effectively in the same way as 
in an ordinary DP process. During the avalanche, several ^4's transmute into -B's. Therefore, once the 
avalanche has stopped, several i?-particles are left behind. These £?-particles in turn perform a simple 
random walk, representing solitary particles in the original model. When two B-particles meet, they may 
trigger a new avalanche of A-particles. The corresponding reaction-diffusion scheme reads 

(1) A^2A, A~>0/B, 2B^A. 

Thus the model consists of two subsystems A and B. Subsystem A is a DP process A «-> 2 A, A — » 
which is coupled via transmutation A — > B to subsystem B. In this subsystem the -B-particles diffuse 
until they annihilate and release a new avalanche of A-particles. Thus, the reaction-diffusion model (^) 
can be interpreted as a cyclically coupled sequence of a DP and a pair annihilation process, as sketched 
in Fig. [| It should be pointed out that Figs. |||| are over-simplified in the sense that they suggest the 
existence of strictly separated dynamic modes. In reality, however, the two modes are not completely 
separated, rather they are entangled and sustain each other in a intricate manner. Nevertheless it is the 
hope of the present study that the AF-process and the simplified reaction-diffusion scheme (Q) exhibit at 
least qualitatively a similar type of nonequilibrium critical behavior. 



4. A LATTICE MODEL FOR CYCLICALLY COUPLED DIRECTED PERCOLATION AND PAIR ANNIHILATION 



In order to study the process (gj) quantitatively, we introduce a three-state model on a square lattice 
with random-sequential updates which is defined by the following dynamic rules: 



reproduction: 
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In the following the numerical analysis will be restricted to the case r = D = 1. Thus, the model is 
controlled by two parameters, namely the rate for particle reproduction p and the transmutation rate r. 
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FIGURE 4. Phase diagram of the lattice model defined in Sec. 
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Figure 5. Simulations at the critical point r = 0.5, p c = 0.6920 starting with random 
initial conditions. Left: The ratio of the particle densities pb/pa tends to a constant. 
Right: Temporal decay of pa and ps multiplied by i 0,21 . 



Obviously, the model has two absorbing states, i.e., the empty lattice and the state with a single diffusing 
S-particlc. The phase diagram in Fig. |j comprises two phases. For low values of p and r, the system is in 
the inactive phase where it approaches one of the two absorbing states. If p and r are sufficiently large, 
an active steady state with non-vanishing particle densities pa and pb exists on the infinite lattice. Here 
we are interested in the critical behavior at the phase transition line. It should be noted that the case 
p = is special. In this case the spreading process does not generate B particles. Hence, starting with 
A-particles, the critical behavior belongs to the DP universality class. For p > 0, however, we observe 
non-DP critical behavior. In fact, an ordinary DP transition seems to be unlikely since the inactive phase 
is charactrized by an algebraic decay of the particle density. Similarly, we can rule out the possibility of 
a DP2 transition since there is neither a ^-symmetry nor a parity conservation law in the model. 
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Figure 6. Left: Finite size data collapse according to the scaling form (g). Right: 
Data collapse for off-critical simulations according to the scaling form (Q) . 

5. Numerical results 



In order to estimate the critical exponents in 1+1 dimensions, we perform Monte Carlo simulations. 
Starting with randomly distributed A's and _B's, we first measure the densities p A {t) and PB(t) up to 3T0 5 
time steps. As can be seen in the left panel of Fig. [s], the ratio ps(i) / ' p A (t) approaches a constant value. 
Assuming an algebraic decay, both quantities should therefore scale with the same critical exponent 

(2) PA (t) ~ p B {t) ~ t~ s . 

The temporal decay of the particle densities is shown in the right panel of Fig. ||. As in the case of 
the AF model, we observe strong corrections to scaling, leading to a considerable curvature of the data. 
However, compared to the AF model these corrections are less severe. Moreover, the curvatures for pa 
and pb have opposite signs. Thus, the local slopes approach the postulated 'true' value of 8 from both 
sides. Seeking for the best compromise, we are able to estimate the critical points by 



T 





0.1 


0.5 


1 


Pc 


0.7674(3) 


0.757(2) 


0.6920(1) 


0.540(1) 



For r = 0.5 the corresponding exponent is given by S = = 0.21(2). Similar measurements for t = 0.1 
and t = 1 (not shown here) suggest that this value is the same for all r > 0. In order to obtain the 
dynamic exponent z — m /u±, we perform finite size simulations at criticality. Here the particle densities 
should obey the scaling form 

(3) PA (t)~ p B (t)~r s f(t/L z ), 

where / is a universal scaling function. Thus, plotting pt s against t/L z , all data sets should collapse onto 
a single curve. Comparing different data collapses it turns out that ps shows a much cleaner scaling 
behavior than p A . As shown in Fig. the best collapse is obtained for i5 = 0.215(15) and z = 1.75(5). 

To obtain the third exponent, we study the behavior pB(t) below and above criticality. According to 
the usual scaling theory for absorbing-state transition, we expect the scaling form 

(4) p B (t)~t- s g(te^), 
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Figure 7. Simulations starting with a single A-particle at the origin. Left: Survival 
probability P(t) and mean square spreading from the origin R 2 (t) averaged over the 
surviving runs. Right: Average number of particles A^i(t) and iVg(t) and their quotient. 



where e = \p — p c \ denotes the distance from criticality. Plotting pB(t)t s against te u ^ (see Fig. the 
best data collapse is obtained for 5 — 0.215(20) and v\\ = 1.8(1). 

In order to cross-check these estimates, we perform dynamic simulations starting with a single pair 
of particles located in the center As usual in this type of simulations, we measure the survival 

probability P(t) that the system has not yet reached one of the two absorbing states, the average numbers 
of particles NA(t) and ATs(i), and the mean square spreading of all particles from the origin R 2 {t) averaged 
over the surviving runs. Assuming that A^(t) and iVs(i) scale asymptotically with the same exponent, 
these quantities should obey the power laws 

(5) P(t)~t- 5 ', N A (t)~N B (t)~t\ R 2 (t)~t 2 l z 

with certain dynamical exponents 5' und rj. As shown in Fig. 0, the survival probability shows a clean 
power law over almost five decades. Moreover, the quotient A^b (t)/A^4 (t) quickly tends to a constant 
value. Fitting power laws to the data shown in Fig. |7[ we obtain the estimates 

(6) 8' = 0.15(1), 77 = 0.21(1), 2/^ = 1.16(4). 

Together with the previous results these estimates satisfy the generalized hyperscaling relation p2{ 

(7) 5 + 5' + rj-d/z = 

within numerical errors. Combining all results, the critical exponents are given by 

p = 0.38(6) , = 0.27(3) , 
( } i/„ = 1.8(1), u ± = 1.0(1). 



It should be noted that the error bars were obtained by assuming power-law behavior. Thus, they do not 
include systematic errors due to possible corrections to scaling emerging after very long time. 
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6. Conclusions 

This work was motivated by recent studies of the annihilation-fission process which exhibits a novel 
type of nonequilibrium critical behavior. Using a scale-invariant space-time plot we have demonstrated 
that critical clusters of the AF process are characterized by an interplay of two different dynamic modes. 
In the high-density mode we observe spreading avalanches whereas the low-density mode is characterized 
by random walks of solitary particles. In order to understand the interplay between these modes from 
a phenomenological point of view, we have introduced an effective model which involves two species of 
particles. The dynamic rules of this model can be regarded as being composed of cyclically coupled 
DP and pair annihilation processes, i.e., it follows the reaction-diffusion scheme A «-> 2 A, A — » B, 
IB — » A. The model exhibits a nonequilibrium phase transition which is in many respects similar to 
the one observed in the AF process. In fact, the AF process and the three-state model proposed in the 
present work have several features in common: 

- They both have two non-symmetric absorbing states, namely the empty lattice and the state with 
a single diffusing particle. 

- There is no unconventional symmetry (such as parity conservation). 

- Both models exhibit a continuous nonequilibrium phase transition with non-DP critical exponents. 

- The visual appearance of critical clusters is very similar (see Fig. [j]). 

Compared to the AF process, the three-state model has several advantages. On the one hand, it is defined 
as a two-site nearest neighbor process. On the other hand, corrections to scaling are less severe, allowing 
us to determine the critical exponents more accurately. 

There are several open questions. Is the critical behavior of the three-state model universal? If this 
is indeed the case, does it represent a independent universality class different from DP and DP2? Are 
the AF process and the present model in the same universality class, i.e., do they have the same critical 
exponents? From the field-theoretic point of view, there is no reason for them to coincide. Yet the 
two classes may 'intersect' in 1+1 dimensions. Very recently, Odor carried out a systematic study of 
the annihilation/fission process performing simulations on a parallel computer combined with generalized 
mean field approximations and coherent anomaly extrapolations [ ^of . He reports two different universality 
classes for low and high values of the diffusion constant d. The corresponding critical exponents are shown 
in Table 1. According to Odor, the exponents for d > 0.5 are in fair agreement with the exponents observed 
in the three-state model. However, this conclusion is still speculative and needs to be substantiated by 
further investigations. 
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